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Abstract 

This paper aims to present the pure field part of the newly developed nonlinear Ex- 
tended Electrodynamics [l]-[3] in non-relativistic terms, i.e. in terms of the electric and 
magnetic vector fields (E,B), and to give explicitly those (3+l)-soliton solutions of the 
new equations which have the integral properties of photons. The set of solutions to the 
new equations contains all solutions to Maxwell's equations as a subclass, as well as, new 
solutions, called nonlinear. The important characteristics scale factor, amplitude func- 
tion, and phase function of a nonlinear solution are defined in a coordinate free way and 
effectively used. The nonlinear solutions are identified through the non-zero values of 
two appropriately defined vector fields T and A4, as well as, through the finite values of 
the corresponding scale factors. The intrinsic angular momentum (spin) is also defined. 
A limited superposition principle (interference of nonlinear solutions), yielding the well 
known classical coherence conditions, is found to exist. 
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1 Introduction 



The 19th century physics, due mainly to Faraday and Maxwell, created the theoretical concept 
of electromagnetic field, i.e. extended (in fact, infinite) object, having dynamical structure. 
The concepts of flux of a vector field through a 2-dimensional surface and circulation of a 
vector field along a closed curve were coined and used extensively. Maxwell's equations in their 
integral form establish where the time-changes of the fluxes of the electric and magnetic fields 
go to, or come from, in both cases of a closed 2-surface and a 2-surface with a boundary. We 
note that these fluxes are specific to the continuous character of the physical object under 
consideration and it is important also to note that Maxwell's field equations have not the sense 
of direct energy-momentum balance relations as the Newton's law p = F has. Nevertheless, 
they are consistent with energy-momentum conservation, as is well known, the corresponding 
local energy-momentum quantities are quadratic functions of the electric and magnetic vectors. 

Although very useful for considerations in finite regions with boundary conditions, the pure 
field Maxwell's equations have time-dependent vacuum solutions (in the whole space) that give 
inadequate models of the real fields. As a rule, if these solutions are time-stable, they occupy 
the whole 3-space or an infinite subregion of it, and they do not go to zero at infinity, hence, 
they carry infinite energy and momentum. As an example we recall the plane wave solution, 
given by the electric and magnetic fields of the form (in a specally chosen coordinate system) 

E = {u(ct + ez),p(ct + ez), oj; B = {ep(ct + ez), — eu(ct + ez), oj, e — ±1, 

where u and p are arbitrary differentiable functions. Even if u and p are soliton-like with respect 
to the coordinate z, they do not depend on the other two spatial coordinates (x, y). Hence, the 
solution occupies the whole TZ 3 , or its infinite subregion, and clearly it carries infinite integral 
energy (we use Gauss units) 

1 f E 2 + B^ 1 f 

W = — / dxdydz — — (u 2 + p 2 )dxdydz = oo. 

An Jn* 2 y 4tt J-r? k ' 

In particular, the popular harmonic plane wave 

u = U cos(cut ± k z .z), p = P sm(ut ± k z .z), c 2 k 2 z = u 2 , U Q = const, P Q = const, 
clearly occupies the whole 3-space, carries infinite energy 

W = — [ (U + P D )dxdydz = oo 

and, therefore, could hardly be an adequate model of a really created field. 

On the other hand, according to Cauchy's theorem for the wave equation OU = (which 
is necessarily satisfied by the components of E and B in the pure field case), every finite (and 
smooth enough) initial field configuration is strongly time-unstable [4]: the initial condition 
blows up radially and goes to infinity (see the next section). Hence, Maxwell's equations 
cannot describe finite and time-stable localized fields. The contradictions between theory and 
experiment that became clear at the end of the last century were a challenge to theoretical 
physics. Planck and Einstein created the notion of elementary field quanta, named later by 
Lewis [5] the photon. The concept of photon proved to be very seminal and has been widely 
used in the 20th century physics. However, even now, after almost a century, we still do not 
have a complete and satisfactory self-consistent theory of single photons. It worths recalling at 
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this place Einstein's remarks of unsatisfaction concerning the linear character of the vacuum 
Maxwell theory which makes it not able to describe the microstructure of radiation [6]. Along 
this line we may also note here some other results and oppinions [7] -[8]. 

In this paper we consider single photons as (3+l)-dimensional extended finite wave-like 
objects, or (3+1) solitary waves, moving as a whole in a consistent translational-rotational 
manner with the speed of light. Their integral characteristics like frequency, period and spin 
are considered to be intrinsically related to their periodic dynamical structure, and the most 
important integral quantity characterizing the single photon seems to be its intrinsic elementary 
action, this being put equal to Planck's constant h. The usually made in quantum theory 
assumption for a point-like character of photons we consider as inadequate, it can never give 
a satisfactory explanation of where the characteristic frequency of a free photon comes from. 
Indeed, the point approximation means structurelessness, so, the frequency may appear only if 
an external force acts on the point particle, but this means that the particle is not free, hence, its 
velocity must not be constant and its spin momentum is not an intrinsic characteristic, which 
contradicts the usual understanding of Planck's formula E = his. In other words, all objects 
that obey the Planck's formula E = hv do not admit point approximation, hence, the solitary 
wave view, considered as a first step, seems much more natural. Extended Electrodynamics 
(EED) [l]-[3] was built mainly to meet this solitary wave view, i.e. to give a consistent field 
description of single photons, and this paper gives a completely non-relativistic approach to 
the pure field part of EED. This gives new insights into the problem. 

Our assumption that single photons are objects of a finite solitary wave nature has as its 
mathematical representation the soliton concept. In accordance with this concept the com- 
ponents of the corresponding E and B at every moment t have to be smooth nonsingular 
functions and different from zero only inside a finite 3-dimensional region Q t C TZ 3 . Moreover, 
a t-periodic process of constant frequency has to accompany the photon's translational motion 
as a whole. Hence, we have to be ready to meet all difficulties coming from the unavoidable 
requirements for using nonlinear partial differential equations having finite with respect to the 
spatial variables (x,y,z) solutions. This goes along with the Einstein's view that "the whole 
theory must be based on partial differential equations and their singularity-free solutions" [9]. 



2 The new equations 

We recall first Maxwell's equations in the pure field case: 

rotB - — = 0, divE = 0, (1) 

rotE + — = 0, divB = 0, (2) 

where £ denotes the product ct, c is the velocity of light in vacuum and t is the time variable. 
From these equations we get the well known Poynting relation 

E 2 + B 2 l mExBy (3) 



d£ 8tt 4tt 

We explain now why Maxwell's equations (l)-(2) have no (3+1) soliton-like solutions. As 
we know from textbooks on Classical Electrodynamics (CED), (e.g. [10]) from (1) and (2) it 
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follows that every component U of E and B necessarily satisfies the wave equation 



UU = U tt - c 2 [U xx + U yy + U zz ] = 0. (4) 

We are interested in the behavior of U at t > 0, if at t — the function U satisfies the initial 
conditions 

dU 

U\ t= Q = f(x, y,z), — = F(x, y, z), 

at t=o 

where the functions f(x, y, z) and F(x, y, z) are finite, i.e. they are different from zero in some 
finite region Q C TZ 3 and have no singularities. Besides, we assume also that / is continuously 
different iable up to third order, and F is continuously different iable up to the second order. 
Under these conditions Poisson proved (about 1818) that a unique solution U(x,y, z;t) of the 
wave equation is defined, and it is expressed by the initial conditions / and F through the 
following formula (a good explanation is given in [4]): 



U(x ' v ' z ' t) = -Lm 



Set 



- f ^-da r ) , (5) 
Js ct r 1 



where P is a point on the sphere S of radius r = ct, centered at the point (x, y, z), and da r is 
the surface element on S r=ct . 

The above formula (5) shows the following. In order to get the solution at the moment t > 
at the point (x,y,z), being at an arbitrary position with respect to the region Q Q , where the 
initial condition, defined by the two functions / and F, is concentrated, we have to integrate 
f/r and F/r over a sphere S r=c t, centered at (x, y, z) and having a radius of r = ct, and then to 
form the expression (5). Clearly, the solution will be different from zero at the moment t > 
only if the sphere S r=ct crosses the region Q Q at this moment. Consequently, if r\ = ct\ is the 
shortest distance from (x,y,z) to Q Q , and r 2 = ct 2 is the longest distance from (x,y,z) to Q Q , 
then the solution at (x,y,z) will be different from zero only inside the interval £2) ■ 

From the point of view of an external observer this means the following. The initially 
concentrated in the region Q Q perturbation begins to expand, it comes to an arbitrary point 
(x,y,z) at the moment t 1 > 0, makes it "vibrate" ( i.e. our devices show the availability of 
field at this point) during the time interval At = t 2 — ti, after this the disturbed point goes 
back to its initial state and our devices find no more field there. Through every point out of Q Q 
there will pass a wave, and its forefront reaches the point (x, y, z) at some moment t\ while its 
backfront leaves the same point at the moment ti> t\. Roughly speaking, the initial condition 
"blows up radially" and goes to infinity with the velocity of light. 

This rigorous mathematical result shows that every finite initial condition induces strongly 
time-unstable free finite time-dependant solution of Maxwell's equations in vacuum, so these 
equations have no finite and smooth enough, i.e. nonsingular, time-dependent solutions, which 
could be used as models of real photons, as viewed by us. 

Hence, if we want to describe 3-dimensional time-dependent soliton-like electromagnetic 
formations (or configurations) it is necessary to leave off Maxwell's equations and to look for 
new equations for E and B. 

On the other hand we know that Maxwell's theory is widely used in almost all natural 
sciences and electrical engineering, so, it does not seem reasonable to leave it off entirely and 
to look for a completely new theory. Moreover, in all energy computations for finite volumes, 
it gives very good results. This suggests, at this stage, to look for some extension of the theory, 
i.e. to extend in a nonlinear way the equations, keeping all solutions to Maxwell's equations 
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as solutions to the new equations and keeping the energy-momentum relations of Maxwell's 
theory as relations of the new theory. In the same time we must incorporate new solutions with 
corresponding to our purpose properties. 

The road we are going to follow in searching for the appropriate nonlinearization of (l)-(2) 
is suggested mainly by two ideas: the idea of local energy-momentum conservation, and the 
idea of invariance of Maxwell's equations (l)-(2) with respect to the transformation 

(E,B)-(-B,E). (6) 

We begin by recalling the second Newton's law in mechanics: p=F, its true sense is local 
momentum balance, i.e. the momentum gained by the particle is lost by the external field. So, 
if the field is absent: F = 0, then the particle will not lose, or gain, momentum, and we get the 
evolution equation for a free particle: p = mv = 0. On the other hand, if the particle is absent 
then the field will not change its momentum and we get the equation F = 0. Usually, F depends 
on the characteristics of the particle, as well as, on the characteristics of the field. In order 
to make possible the interpretation of the relation F = as a pure-field equation we have to 
express, if possible, this F in terms of the corresponding field functions (and their derivatives) 
only. Fortunately, this can be done in electrodynamics, Maxwell's eqiations in presence of 
sources make it possible. In this way we shall obtain explicitly in terms of the field functions 
and their derivatives one of the possible expressions describing how much momentum the field 
is potentially able to transfer locally to another physical system in case of its presence at the 
corresponding point and its ability to absorb this momentum, i.e. to interact with the field. 
Of course, the field may have various such momentum (and corresponding energy) exchanging 
abilities. We'd like to note at this moment that the following considerations (up to equation 
(7)) have a suggestive nature only, they do not prove what we are going to assume finally as 
the new equations. 

In order to carry out the above idea we recall first that if the other system consists of 
charged particles, as it is in Maxwell theory, the corresponding force is the well known Lorentz' 
force, acting on a particle of electric charge e: F = eE + ^(v x B). In case of a continuous 
distribution of particles with charge density of p and current j = pv the Lorentz' force is 

F = pE + -(v x B) = pE + -(j x B). 

c c 

The corresponding Maxwell's equations with non-zero charge distribution p and current j = pv 
in this case read: 

„ dE 4tt. , „ A 
rotti — — = — j, divE = 4np. 
d£ c 

These last equations make possible to substitute j and p into the above given Lorentz' force. 
Having this done we put F = and then we forget about the character of the "other system" 
(charged particles). We interpret now F, so obtained and expressed through the field functions 
and their derivatives only, as a definite quantity of momentum which the field is potentially able 
to transfer (locally) to any (continuous) physical system that is able to absorb it. This defines 
quantitatively one of the field's momentum exchanging abilities. 

The above suggestive considerations make us assume our first extended equation: 

x B + EdivE = 0. (7) 
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This vector equation (7) extends Maxwell's pure field equations (1) in the sense that (7) implies 
no more (1), i.e. (7) may have solutions which do not satisfy (1) and these new solutions, at 
least some of them, are considered as admissible from physical point of view, i.e. describing 
some physical reality. The physical sense of (7) is quite clear: no field momentum is lost in this 
way. At the same time relation (7) describes some internal, i.e. between E and B, redistribution 
of the field energy-momentum during the field's time evolution. The nonlinearity of (7) is also 
obvious. We'd like to stress once again the suggestive character of the obove considerations, 
equation (7) is an assumption, we do not consider it as a consequence of Maxwell equations 
with non-zero current, on the contrary, we consider it as a pure field equation. 

Now we look for another momentum exchanging ability of the field, different from the one 
given by the left-hand side of (7), and directed in general to new physical systems. In order to 
come to such an expression we make use of the above mentioned (E, B) — > (— B, E) invariance 
of the pure field Maxwell's equations, known as electro-magnetic duality. This invariance is 
valid also for the energy w and momentum S (i.e. Poynting's vector) densities: 

E 2 + B 2 B 2 + E 2 - c c c 
w = ► , S=-ExB-4 BxE=— ExB. 

8tt 8tt 4tt 4tt Att 

Moreover, the basic energy- momentum balance relation of Poynting (3), is also invariant with 
respect to this transformation (6). This suggests that, transforming (7), i.e. replacing in (7) E 
by — B and B by E, we should obtain also a true and valuable relation, since (7) describes now 
intra-field local momentum balance. In this way we obtain our second vector equation: 

(<9B \ 
rotE + —j x E + BdivB = 0. (8) 

The left-hand side of (8) defines explicitly another momentum exchanging ability of the field, 
and relation (8) defines another way of internal field energy-momentum redistribution with 
time. Note that (8) is obtained from (7) in the same way as (2) is obtained from (1). 

We complete this process of extension of Maxwell's equations by adding two new invariant 
with respect to the same transformation (6) equations, which also have the physical sense of 
intra-field local energy-momentum balance: 

rotB - — \ x E - EdivB - BdivE = 0, (9) 

B (rotB - " E. (rotE + ^ = 0. (10) 

Remark: In the relativistic formulation of EED [3] the two relations (9) and (10) have a natural 
interpretation of energy-momentum transfers between F and *F, where F is the conventional 
electromagnetic field tensor constructed by means of E and B, and *F is its dual (constructed 
similarly by means of — B and E). This interpretation says that these transfers are mutual: 
F # *F, and always in equal quantities. In other words, the relativistic formalism considers 
the couples (E, B) and (— B,E), or F and *F, as two componenents of a new more general 
mathematical object [2]- [3]. 

Note that under the transformation (6) equations (7) and (8) transform into each other, 
while equations (9) and (10) are kept the same (up to a sign of the left-hand side). 
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Equations (7)-(10) constitute our new system of equations for the electromagnetic field in 
vacuum. Obviously, they do not introduce new parameters, they are non-linear and all solutions 
to Maxwell's vacuum equations (l)-(2) are solutions to the new equations (7)-(10). We are going 
now to study those new solutions to (7)-(10), which satisfy the conditions 

<9B <9E 
rotE + — ^ 0, rotB - — ^ 0, divE ^ 0, divB ^ 0. (11) 
of of 

For further convenience, all solutions to (7)-(10), which satisfy (11), will be called nonlinear. 



3 Properties of the Nonlinear Solutions 

The first two, almost obvious, properties of the nonlinear solutions follow directly from (7) and 
(8) and are given by the relations 

E.B = 0, (12) 

(<9B \ / dE \ 

rotE + — J = 0, E. (rotB - — J = 0. (13) 

Relation (12) says that the electric and magnetic vectors of every nonlinear solution are nec- 
essarily orthogonal to each other at every point, the algebraic property that Maxwell's (linear) 
plane wave solution has. From (13) it follows that the Poynting's relation (3) is true for all 
nonlinear solutions, and this justifies the usage of all energy- momentum quantities and relations 
from Maxwell's theory in the set of nonlinear solutions of the new equations. We can consider 
the left-hand sides of relations (13) as the energy quantities which the field is potentially able 
to transfer to some other physical object. We note also the obvious invariance of (12)-(13) with 
respect to transformation (6). 

We are going to show now that all nonlinear solutions satisfy the relation 

E 2 = B 2 . (14) 

In order to prove (14) let's take the scalar product of equation (8) from the left by B. We 
obtain 

B.j(VotE + x e| +B 2 divB = 0. (*) 
Now, multiplying (9) from the left by E and having in view (12), we obtain 

E.j (rotE + x B j - E 2 divB = 0. 
This last relation is equivalent to 



-B.j (rotE + x E j - E 2 divB = 0. 



Now, summing up (*) and (**), in view of divB 7^ 0, we obtain (14). 

Relation (14) is also true for the linear plane electromagnetic wave. It requires for all 
nonlinear solutions a permanent equality of the energy densities carried by the electric and 



7 



magnetic fields, although permanent mutual energy- momentum flows run between E and B, 
which means that these two flows are always in equal quantities. 
We introduce now the following two vector fields: 

- _ <9B E x B 

T = rotE + — + ^-^divB, (15) 
<9E E x B 

M = rotB - — - — — —divE. (16) 
<9£ |E x B| 

It is obvious that on the solutions of Maxwell's equations (l)-(2) T and M are equal to 

— * — * — * — * 

zero. Note also that under the transformation (6) we get T — > — M. and A4 — > T . We shall 
show now that on the non-zero nonlinear solutions of our equations (7)-(10) T and M. are 
colinear to E and B respectively. Indeed, consider the products T x E and M. x B. Since 

(E x B) x E = -E x (E x B) = — [E(E.B) - B(E.E)] = B|E| 2 

and |E x B| = E 2 = B 2 , we obtain (see eqn.(8)) 

/ dB\ 
fxE= rotE + — x E + BdivB = 0. 

— * 

In the same way we get M. x B = 0. Hence, we can write the relations 

f = f 1 .E, M = f 2 .B, (17) 

where f\ and f 2 are two functions, and further we consider the interesting cases /i ^ 0, oo; 
f 2 7^ 0, oo. Note that the physical dimension of fi and f 2 is the reciprocal to the dimension of 
coordinates, i.e. [fi] = [f 2 ] = [length] -1 . Note also, that JF and M. are mutually orthogonal. 

We shall prove now that f\ = f 2 . In fact, making use of the same formula for the double 
vector product, used above, we easily obtain (see eqn.(9)) 

f xB + MxE = 

( dB\ ( dE\ 

= IrotE + — J x B + IrotB - —j x E - EdivB - BdivE = 0. 

Therefore, 

fxB+MxE= 
= fiE x B + f 2 B x E = (/i - f 2 )E x B = 0. 

— * — * 

The assertion follows. Now from (14) and (17) it follows also the relation \F\ = \M\. 
Definition 1. The quantity 

1 1 |E| |B| , x 

L ( E ' B ) = ITT = ITT = wi=nrr ^ 

\h\ \h\ {J 7 ] \M\ 

will be called the scale factor for the nonlinear solution (E,B). 

Obviously, L(E,B) = L(— B,E), and for all non-zero nonlinear solutions we have 
< L < oo, while for the linear solutions L — > oo. 
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4 Photon-Like Solutions 



As we mentioned earlier, we consider photons as finite nonsingular objects moving translationaly 
along straight lines in the 3-space with the velocity of light. The direction of motion is assumed 
to be that of (E x B). This means that the integral curves of (E x B) have to be straight lines. 
If we choose the coordinate system (x, y, z) so that this direction of the translational motion to 
coincide with the direction of the coordinate line z, in this coordinate system the vector fields 
E and B will have non-zero components only along x and y: 

E\ = (u,p, 0); B = (m, n, 0), 

so, (E x B) = (0, 0, un—pm). Now from E.B = and E 2 = B 2 it follows m = ep and n = —eu, 
e — ±1. Hence, 

E = (u,p,0); B = (ep, -eu, 0), E x B = [0, 0, -e{u 2 + p 2 )], 

and we have to determine just the two functions u and p. 

Let's substitute these E and B into the left hand sides of equations (7)-(10). We obtain: 

rotE + — J x E + BdivB = [0, 0, ep(p(. - ep z ) + eu(u^ - eu z )\; 

(dB \ 
rotB — — J x B + BdivB = [0, 0, eu(u^ — eu z ) + ep(p^ — ep z )\; 

(dB \ 
rotB + — J x B = [eu(p x - u y ),ep(p x - u y ), -ufa - ep z ) + p{u^ - eu z )\; 

{rotB - x E = [ep(u x + p y ), eu(u x + p y ), -p(u^ - eu z ) + - ep z )\; 

—BdivB = [-eu(p x - u y ), -ep(p x - u y ), 0]; 
-BdivB = [-ep(u x +py),eu(u x +p y ),0\; 

(dB \ 

rotB - — j = -ep(u£ - eu z ) + eu(p^ - ep z ), 

where the indices of u and p mean the corresponding derivatives. 

It is seen that equations (9)-(10) are satisfied identically, and equations (7)-(8) reduce to 
only one equation, namely 

u( Ui . - eu z ) + p(p£ - ep z ) = i[(w 2 + p 2 )^ - e{u 2 + p 2 ) z ] = 0. (19) 
The solution to this equation is 

u 2 +p 2 = <p 2 (x,y,i + ez), (20) 
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where is an arbitrary differentiable function of its arguments. This relation means that the 
energy density 

JL(E> + B>) = -L( U ' + tf 

is a running wave along the coordinate z. Hence, our equations determine the field components 
u and p up to a bounded function <p(x,y, z,£), \(p\ < 1: 



u 



Reversely, 



'y.tp, p = ±<i>.y/l - ip 2 . 
= ±\Ju 2 + p 2 , ifi 



u 



(21) 



(22) 



±\ / u 2 + p' 2 

The above relations show that instead of u and p we can work with and ip. Equations (7)- 
(10) require only to be running wave along z (in this coordinate system), and ip to be bounded 
function. In all other respects these two smooth functions are arbitrary. Hence, they may be 
chosen finite with respect to the spatial coordinates (x,y,z). Hence, the nonlinear equations 
(7)-(10) allow (3 + 1) soliton-like solutions. Note that, since ip is bounded, it is sufficient to 
choose just to be spatially finite. 

We are going to show now that the two functions and ip have a certain invariant sense 
and can be introduced in a coordinate free way. 

First, let's denote by a the invariant 



a 



Since E.B = and E 2 = B 2 we have a 



IEI 



|B| 1 , and the local frame 



X = («E, -aeB, -a 2 eE x b) 



is orthonormal. This frame is defined at every point where the field is different from zero. At 
every point we have also the frame of unit vectors generated by the coordiate system chosen: 



| = (1,0,0), e^| = (0,l,0), e. = | 



(0,0,1) 



We represent now the frame vectors of x through the frame vectors of x°, and obtain the matrix 



A 



\J V? +p' 2 yj u 2 +p 2 





yfU 2 +p 2 y/U 2 +p 2 

1 



This matrix has three invariants: l\ = trA, li—tht sum of all principal minors of second order, 
h = detA. We find 



2u 



\fu 2 +p 



= + 1; I 3 = detA = 1. 



Clearly, \{h - 1) < 1. 



Hence, we can define and ip through a and these invariants in the following way: 

<P = ± x Ja- 2 I 3 (A), <p=UtrA-\). (23) 
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Definition 2. The functions and ip, defined by (23) will be called amplitude function and 
phase function of the corresponding nonlinear solution, respectively. The function arccos(p) 
will be called phase of the solution. 

— * — * 

For T and M. we obtain 



F = [ £ (Pt ~ £ Pz), -e(«{ - eu z ), 

M = - eu*), -(p { - ep z ), 

We shall express the scale factor L through and ip. We obtain 



\F\ = \M\ 



Therefore, since |E| = |B| = |0|, the scale factor L is obtained as function of p and its first 
derivatives only, 



|E| 



B 



(24) 

\M\ \<Pz~e<Pz\ 

Now we shall separate a subclass of nonlinear solutions, called almost photon-like, through 
the following conditions on ip and L: 



dip dip n dp dip 

M 7r + ^7r = ' Pif- U ir 

dx dy dx dy 



o- d ± 



0. 



(25) 



The invariant sense of the first three equations of (25) is that the phase function p is a first 
integral of the vector fields E, B, E x B. From the third equation of (25) is clearly seen that in 
this coordinate system p does not depend on z. The first two equations of (25), considered as 
an algebraic linear homogeneous system with respect to the two derivatives, yield p x = ip y = 
because the corresponding determinant is always nonzero: u 2 + p 2 ^ 0. Hence, ip may depend 
only on £. In view of (24) and ip z = the fourth equation §| = means L = constant. Hence, 
relation (24) turns into equation for ip: 



VT 



p< 



The obvious solution is (both signs on the right are admissible) : 



ip = COS (kj- + P 



COS (^jt + Po^j , 



(26) 



(27) 



where k = ±1, f3 Q = const. Since ip is a periodic function with respect to t, then c/L has the 
physical interpretation of frequency, and this frequency has nothing to do with the concept of 
frequency in classical electrodynamics since it is defined by L, and L is not defined in Maxwell's 
theory. 

It is clearly seen the consistent translational-rotational behavior of the solution obtained: 
the electric and magnetic vectors 



E = 



cos (k^ + fi^j , (j) sin (k^ + p o ^j , 
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e<j) sin [k^ + (3^j , -e<f) cos (k^ + (3^ , 



run along z: <fi = (f)(x,y,^ ± z), and rotate (left or right: k = ±1) with the frequency 
u = (c/L) = l/T. 

In order to separate the photon-like solutions we recall that the photon's characteristic 
quantity is its integral intrinsic angular momentum, or spin, being equal to the Planck's constant 
h. Namely h represents quantitatively in a unified manner the rotational and translational 
aspects of its dynamical nature: for all photons the product WT has the same value h, although 
W and T may be different for the different photons. That's why every photon should be able 
to determine its own scale factor L = const in order to have a cosine periodic phase function 
and to obey the Planck's law: h = WT = WL/c. The photon's intrinsic periodic process 
demonstrates itself in our approach through the (left or right) rotation of the pair (E,B). 
Since these two vectors are orthogonal to each other and with equal modules: |E| = |B|, the 
basic local quantity appears to be the area of the square defined by E and B at every point, and 
this area is equal to |E x B|. During one period T this square performs a full rotation around 
the direction of propagation and this gives the local action |E x B|.T. In order to obtain the 
integral T-action of the solution we have to sum up all these local actions. 

The above described idea is easily represented mathematically. In fact, let u, p and u 2 + p 2 
be spatially finite functions. Then the integral energy 

W = — [ (u 2 + p 2 )dxdydz < oo 

All J-r 3 

is finite. For every almost photon-like solution we define the local spin vector S by 

S = = llli! L = const ^ o < L < oo. (28) 

Now, the integral intrinsic action, or integral spin S of the solution, is defined by 

S = J \S\dxdydzdt. (29) 

We obtain 

S = WT. (30) 

We note once again that this approach works because W = const < oo, < L = const < oo 
and the solution is soliton-like, i.e. it is finite, it has periodic dynamical structure and is time- 
stable. Clearly, no solution of Maxwell's equations (l)-(2) in the whole space has all these 
properties. 

Definition 3. A nonlinear solution will be called photon-like if it is spatially finite, if it 
satisfies conditions (25) and if its integral spin S is equal to the Planck constant h: S = h. 

Finally, we consider briefly the problem of interference of photon-like solutions: if we have 
two photon-like solutions 



Ei = 



(fr 1 (x,y,£ + e 1 z)cos +Pij ,<j>i(x,y,£ + eiz)sm + Pij ,0 
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B 1 



£i0i(x, y, £ + £i2:)sin J- + ft^ - £i0i(x, y, £ + £i2:)cos ^1 J- + Pij , 

i 2 (x, y, £ + £ 2 z)cos ^2-^- + P^j , 02 (a:, y, £ + £ 2 ^)sin ^k 2 J- + , 
£ 2 2 (x,y,£ + £ 2 z)sin («2-|- + /? 2 j - £ 2 2 (:r, y, £ + £ 2 ^)cos f^J- + /3 2 j ,0 



we ask: under what conditions their sum (Ei + E 2 , B x + B 2 ) will be again a nonlinear solution? 
Having done the corresponding elementary computations, we come to the following important 
conclusion: if 

£i = £2, «i = «2, -^1 = ^2 (31) 

then the sum (Ei + E 2 , Bi + B 2 ) is again a nonlinear solution. 

These relations (31) coincide with the well known from CED coherence conditions. They 
say that the two photon-like solutions will interfere, i.e. their sum will be again a solution, if: 

1. They propagate along the same direction: E\ — £ 2 , 

2. They have the same polarization: k\ — « 2 , 

3. They have the same frequency: v\ = i/ 2 , i.e. the same scale factors. 

Recall that in CED these coherence conditions do not follow directly from the theory as 
necessary conditions, e.g., condition 3 requires some time averaging. Moreover, CED is a linear 
theory and the sum of any two or more solutions is again a solution, for example, the sum 
"plane wave + a spherically symmetric field" is again a solution but no interference features 
are available. 

In EED, which is a non-linear theory and there is no superposition principle in general, the 
interference of photon-like solutions is a very special case and it is a remarkable result that 
the experimentally found coherence conditions (31) appear from the nonlinear equations as 
necessary conditions, otherwise the sum will not be a solution. 

The computation shows that the energy density w of the sum-solution is given by 



w 



47T 



(0i) 2 + (0 2 ) 2 + 2</. 1 (/. 2 cos(/5 1 -/5 2 ) 



(32) 



and this relation (32) allows to talk about interference instead of superposition: in our approach 
the allowed superposition of photon-like solutions leads always to interference. 

Of course, from the non-linear point of view, these interference phenomena are of some 
interest only for those soliton-like solutions which at a given moment occupy intersecting re- 
gions, otherwise the interference term 20i0 2 cos(/3i — /3 2 ) in (32) is equal to zero. Since the 
two summonds follow the same direction of motion as a whole, then the sum-solution will be a 
time-stable solution, but it will not be photon-like one. 



5 Conclusion 

This paper presents a non-relativistic formulation of an extension of the pure field Maxwell 
equations. The main purpose of the extension is to give a mathematical description of the 
viewpoint that photons are extended but finite objects and that their existence is based on a 
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joint and consistent translational-rotational internal dynamics: a straight-line (translational) 
motion as a whole with the velocity of light and a rotation of the mutually orthogonal and 
perpendicular to the translation electric and magnetic vectors. The general mathematical 
concept of (3 + l)-solitary wave, or soliton, was found to be the most adequate one to this 
physical notion. Our nonlinear equations (7)-(10) for the couple of vector fields (E, B) realize 
directly the idea for local energy-momentum conservation, i.e. the idea that the free field 
dynamics (or time evolution) is determined by a definite and permanent intra-field local energy- 
momentum redistribution. Compare to Maxwell's approach, this is a new look on the field's 
evolution defining equations, it carries the idea of the second Newton's law in mechanics to 
continuously distributed in the 3-space physical systems (electromagnetic fields). It is also a 
new moment that the field has more than one potential abilities to exchange energy-momentum 
with other physical systems. This nonlinear approach turned out to be successful in view 
of the existence of appropriate photon-like solutions. Every photon-like solution has finite 
integral energy W, has its own scale factor L = const, phase function if of cosine type and 
corresponding to L and <p frequency v — cj L, or period T = L/c. Every photon-like solution 
carries intrinsic angular momentum of integral value equal to the Planck's constant h = WT, 
which is the famous Planck's formula. The equations (7)-(10) and the additional conditions (25) 
do not determine the spatial structure of the solution and it should be so, because it is hardly 
believable that all photons must have the same shape, structure and extension. Moreover, 
if the amplitude function consists of many non-overlaping appropriate bumps, we obtain a 
mathematical image of a number of coherent photons, i.e. a many-bump photon-like solution, 
and the corresponding integral spin momentum will be nh, where n gives the number of bumps. 
Clearly, large enough parts of the classical plane wave can be considerd as macro-pictures of 
such a many-bump photon-like solution. Nonlinear solutions with intrinsic angular momentum 
different from h, or nh, where n = 1,2,..., are also, in principle, allowed. The remarkable 
limited superposition principle permits interference in the frame of photon-like solutions only 
if they are coherent in the sense of (31). 

The existence of localized photon-like solutions in the pure field case suggests to make the 
corresponding extension of Maxwell's equations in presence of external fields (media) and to 
look for (3+l)-localized solutions with non-zero proper mass. Such an extension in relativistic 
terms was made and published [2], where a large family of variously shaped (3 + 1) soliton-like 
solutions with non-zero mass and well defined conserved quantities is also given. 

We would like to acknowledge the substantial support of the Bulgarian National Fund 
"Science Research" through Project — 718 in doing this research. 



14 



REFERENCES 



[1]. S.Donev, M.Tashkova, Proc.R.Soc.of Lond., A 443, (301), 1993. 

[2]. S.Donev, M.Tashkova, Proc.R.Soc. of Lond., A 450, (281), 1995. 

[3]. S.Donev, M.Tashkova, Annales de la Fondation Louis de Broglie, vol.23, No. No. 2, 3 
(1998). 

[4]. S.J.Farlow, "Partial Differential Equations for Scientists and Engineers", John Wiley 
& Sons, Inc., 1982. 

[5]. G.N.Lewis, Nature 118, 874, 1926. 

[6]. A.Einstein, Sobranie Nauchnih Trudov, vols. 2, 3, Nauka, Moskva, 1966. 
[7]. M.Planck, J.Franklin Institute, 1927 (July), p. 13. 

[8]. J.J.Thomson, Philos.Mag.Ser. 6, 48, 737 (1924), and 50, 1181 (1925), and Nature, 
vol.137, 23 (1936); N.Rashevsky, Philos.Mag. Ser.7, 4, 459 (1927); W.Honig, Found.Phys. 
4, 367 (1974); G. Hunter, R.Wadlinger, Phys.Essays, vol.2,158 (1989). 

[9]. A.Einstein, J.Franklin Institute, 221 (349-382), 1936. 

[10]. J.D.Jackson, CLASSICAL ELECTRODYNAMICS, John Wiley and Sons, Inc., New 
York-London, 1962. 



15 



